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Abstract 

A search for excited fermions f * of the first generation in e^p scattering at the col- 
lider HERA is presented using HI data with an integrated luminosity of 37 pb^^. All 
electroweak decays of excited fermions, f* —>■ f^,fW,iZ are considered and all possible 
final states resulting from the Z or W hadronic decays or decays into leptons of the first 
two generations are taken into account. No evidence for f * production is found. Mass 
dependent exclusion limits on cross-sections and on the ratio of coupling constants to the 
compositeness scale are derived. 
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1 Introduction 



Models of composite leptons and quarks [|I]] were introduced in an attempt to provide an ex- 
planation for the family structure of the known fermions and for their pattern of masses. A 
natural consequence of these models is the existence of excited states of leptons and quarks. It 
is often assumed that the compositeness scale might be in the TeV region, which would give 
excited fermion masses in the same energy domain. However, the dynamics at the constituent 
level being unknown, the lowest excitation states could possibly have masses of the order of 
a few hundred GeV. ElectronQ-proton interactions at very high energies provide an excellent 
environment to look for excited fermions of the first generation. 

In this paper a search for excited fermions is presented using e^p HERA collider data of the 
HI experiment. The data collected from 1994 to 1997 at positron and proton beam energies of 
27.5 GeV and 820 GeV respectively correspond to an integrated luminosity of 37 pb^^. The ex- 
cited fermions are searched for through all their electroweak decays into a fermion and a gauge 
boson. The subsequent W and Z gauge boson decays considered are those involving electrons, 
muons, neutrinos or jets. This analysis profits from an increase in statistics by more than a 
factor of 10 compared to previous HI searches ^, and more than a factor of 4 compared to 
published results by the ZEUS collaboration . 

The paper is organized as follows. The phenomenological model used to interpret the results 
of the search for excited fermions is discussed in section ^. The HI detector and the data pres- 
election criteria are described in section ^. The generators used for the Monte Carlo simulation 
of the Standard Model events and excited fermion signals are briefly presented in section^. The 
analyses for the various possible final state topologies are described in section |^. The search 
results are interpreted in section ^ and a summary is presented in section [7|. 

2 Phenomenology 

Excited electrons (e*) could be singly produced in ep collisions through t— channel 7 and 
Z boson exchange (fig. Single production of excited neutrinos {u*) could result from 
t— channel W boson exchange (fig. p?). In the same way excited quarks (g*) could be produced 
through t— channel gauge boson exchange between the incoming positron and a quark of the 
proton (fig. |T]c). 

Excited fermions can carry different spin and isospin assignments In some models 
quarks and leptons are composites of a scalar and a spin | constituent and the lowest lying exci- 
tations have spin |. Alternatively excited fermions could consist of three spin | constituents and 
in this case the lowest lying excitation levels could appear as spin | states ||]. Given that the 
lowest spin states are often considered as the most probable, we use a model [0, |^, ^ in which 
excited fermions are assumed to have spin | and isospin |. This model describes the interaction 
between the excited fermionic particles f *, the gauge bosons and the ordinary fermionic matter 
f by an effective Lagrangian. Both left-handed and right-handed components of the excited 

^The term 'electron' stands generically for electron or positron. 
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Y, Z,W, g 



Figure 1: Leading diagrams for the production and decay of excited fermions in ep collisions. 



fermions form weak isodoublets F£ and F^. In order to protect the light leptons from radia- 
tively acquiring a large anomalous magnetic moment [|T^ [TT]], couplings of excited fermions 
to ordinary fermions of both chiralities should be avoided. We choose to consider couplings 
to left-handed fermions only, in which case only the right-handed component of the excited 
fermions takes part in this magnetic-moment type Lagrangian. The form of the Lagrangian is: 



int 



2 G;^]Fl + h.c. 



(1) 



where W^^, B^,y and G",^ are the field-strength tensors of the SU{2), U{1) and SU{3)c gauge 
fields, t"', Y and A" are the Pauli matrices, the weak hypercharge operator and the Gell-Mann 
matrices, respectively. The g, g' and gs are the standard electroweak and strong gauge couplings. 
A is the compositeness scale and the /, /' and fs can be viewed as form factors (reduced here to 
parameters) allowing for the composite fermion to have arbitrary coupling strengths associated 
to the three gauge groups. 

In this model the coupling constants of the f * with the boson and the fermion f are related 
to the parameters / and /' as: 



CZi*l 



CWf*i 



1 Y 



■ if h cot 9 
f 



w 



Y 

/'—tan 6*1^) 



(2) 



2^2 sin Ow 

where is the third component of the i so spin of the fermion and 9\y is the Weinberg angle. The 
partial widths for the various electroweak decay channels of an excited fermion f * in a fermion 
f and a real boson V are given as 



r(f* ^ fv) 



Ml 2 



1 - 



Ml 
Ml 



2Ml 



(3) 



where Mf . is the mass of the excited fermion. My the mass of the electroweak boson and a the 
fine structure constant. For the excited quark, the partial width to decay to a quark and a gluon 
is obtained replacing {a) by (4/3 Ug) and (cyf.f) by (1/2 fg), where is the strong coupling 
constant. For Mf* values between 50 to 250 GeV and A = Mf*, the intrinsic widths of the 
excited electrons are typically of the order of some hundred MeV for / and /' values ~ 1 while 
in the q* case for ~ 1 it varies between 1 GeV to ~ 10 GeV. 



5 



3 The HI Detector 



A full description of the HI detector can be found in reference [[121] . Here we describe briefly 
the components relevant to this analysis. The interaction region is surrounded by a system of 
drift and proportional chambers covering the angular range 7° < 9 < 176°. The tracking system 
is placed inside a finely segmented liquid argon (LAr) calorimeter covering the polar angular 
range A° < 9 < 154°. The electromagnetic part is made of lead/argon and the hadronic part 
of stainless steel/argon [[T3|]. Energy resolutions oi ge/ E ~ 12%/\fE © 1% for electrons and 
geIEi rr^i 50%/\/E' © 2% for hadrons have been obtained in test beam measurements [ [T^ , [T3| ]. 
A lead-scintillating fiber calorimeteiQ [jl^ is located in the backward^ region (154° < 9 < 
178°) of the HI detector. The tracking system and calorimeters are surrounded by a super- 
conducting solenoid, producing a uniform magnetic field of 1 . 1 5 T in the z direction, and an iron 
yoke instrumented with streamer tubes. Leakage of hadronic showers outside the calorimeter is 
measured by analogue charge sampling of the streamer tubes with a resolution [jTSj] oi as/ E ~ 
100%/V^. Muon tracks are identified from their hit pattern in the streamer tubes. 



4 Event Generators 

Final states of events selected in this analysis contain either a high energy electron (or photon) 
or jets with high transverse energy (or missing transverse momentum). The main backgrounds 
from Standard Model processes which could mimic such signatures are neutral current deep 
inelastic scattering (NC DIS), charged current deep inelastic scattering (CC DIS), photoproduc- 
tion processes, QED Compton scattering and W and Z production. 

For the determination of the NC DIS contributions we used two Monte Carlo samples with 
different modelling of the QCD radiation: 



The first one was produced with the event generator DJANGO [|19p which includes QED 
first order radiative corrections based on HERACLES [^0]]. QCD radiation is imple- 
mented using ARIADNE [|l]] based on the Colour Dipole Model (CDM) [||]. This 
sample, with an integrated luminosity of more than 10 times the experimental luminosity, 
is chosen for the estimation of the NC DIS contribution unless explicitly stated otherwise. 



The second sample was generated with the program RAPGAP [EBH, where QED first 



order radiative corrections are implemented as described above. RAPGAP includes the 
leading order QCD matrix element and higher order radiative corrections are modelled 
by leading-log parton showers. This sample of events corresponds to an integrated lumi- 
nosity of about 2 times the experimental one. 



For both samples the parton densities in the proton are taken from the MRST p4l ] parametriza- 
tion which includes constraints from DIS measurements at HERA up to squared momentum 

^This detector has replaced in 1995 a conventional lead-scintillator sandwich calorimeter [p^. 
■^The forward direction, z >0, from which the polar angle 9 is measured is the proton beam direction. 
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transfer = 5000 GeV^ p5|, and the hadronisation is performed in the Lund string frag- 



mentation scheme by JETSET pTQ. 



The modelling of the CC DIS process is performed by DJANGO using MRST structure 
functions. QED radiation from quark lines is not fully included in the NC DIS and CC DIS 
simulations. Whilst inelastic wide angle bremsstrahlung (WAB) is treated in the generator 
DJANGO, elastic and quasi-elastic WAB (QED Compton scattering) is simulated with the event 



generator EPCOMPT [ p8| , p9| ]. Direct and resolved photoproduction processes (7p), includ- 
ing prompt photon production are simulated with the PYTHIA [ PD[ ] event generator. Other 
processes, corresponding to much lower cross-sections, such as lepton pair production or W 
production have also been simulated. The lepton pair production (77) is simulated using the 



LPAIR generator [|31|]. It should be noted that this generator contains only the Bethe-Heitler 
77 process. However the number of events with two high Ex electromagnetic clusters given 
by LPAIR and by a generator taking into account all (electroweak) tree level graphs and addi- 
tional first order radiative corrections [ ]32[ ] agrees within 5%. The W simulation is made with 
the EPVEC program The luminosities generated for these Monte Carlo simulations vary 



between 3 times to 100 times the experimental luminosity. 

Monte Carlo simulations of excited fermion production and decay are necessary to evaluate 
acceptance losses due to selection requirements. The excited fermion analyses are based on the 
phenomenology described in section 0. The excited lepton (/*) simulation is performed by the 
COMPOS [0] generator which makes use of the cross-section formulae given in reference 
The excited quark generation is done following the cross-section given in reference [p. In both 
cases initial state radiation of a photon from the incoming electron is generated. The photon 
is taken to be coUinear with the electron, with an energy spectrum given by the Weizsacker 
Williams formula. The hadronisation is performed here also by the Lund string model and 
the MRST parametrization of the parton densities is used. The narrow-width approximation is 
assumed and the production and decay of the excited fermions are assumed to factorize. 

All Monte Carlo generators are interfaced to a full simulation of the HI detector response. 



5 Event Selection and Comparison with Standard Model Ex- 
pectation 

In this section the description of the selection criteria for the analyses of the various decay 
channels is organized according to the experimental signatures of the final states. Other details 
of the analyses can be found in reference [|3^]. 

In common for all analyses, background not related to e^p collisions is rejected by requiring 
that there is a primary vertex within ±35 cm of the nominal vertex value, and that the event 
time, measured with the central tracking chamber, coincides with that of the bunch crossing. In 
addition topological filters against cosmic and halo muons are used. A small number of cosmic 
and halo muons are finally removed by a visual scan. 

The identification of electrons or photons, performed in the LAr calorimeter, first relies 
on calorimetric information by exploiting the shape of the energy density expected from the 
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development of an electromagnetic shower to define electromagnetic clusters. An electron is 
identified as an electromagnetic cluster with a track linked to it. A photon in contrast should 
have no track pointing to it within a distance of 40 cm. In some analyses electrons and pho- 
tons are not distinguished and, in this case, only electromagnetic (em) clusters are required. 
Hadronic jets (denoted jets in the following) are searched for in the laboratory reference frame. 



using a cone algorithm adapted from the LUCELL scheme from the JETSET package [pTQ, with 
a radius R = \J brj^ + 50^ = 1, where ri is the pseudorapidity and the azimuthal angle. A 
muon is identified as a well measured central track linked geometrically to a track in the muon 
system or an energy deposit in the instrumented iron. A muon candidate should also satisfy an 
isolation criterium imposed in the pseudorapidity-azimuth [j] — (p) plane by requiring that the 
distances of the muon track to the nearest hadronic jet and to the closest track be greater than 1 
and 0.5 in i] and (p respectively. 

The event selection makes use of the global variables (E — P^), P^^° and P™'"-*- described 
in the following. 



• {E—Pz) = J2i {Ei—Pz. ) where Ei and P^^ = Ei cos 9i are the energy and the longitudinal 
momentum measured in a calorimeter cell i. For an event where the only particles which 
remain undetected are close to the proton direction, momentum conservation implies that 
{E — Pz) nearly equals twice the energy of the incoming positron (55 GeV). 

• p™^° = I pmio I \y]^ej-e pmio ^ missing transverse momentum vector with components 
calculated by summing over all energy deposits recorded in cells of the LAr and back- 
ward calorimeters. For the study of the channels including a muon in the final state, this 
sum is extended to the energy deposits in the instrumented iron. This P|."'° variable mea- 
sures the transverse energy of undetected particles (neutrinos) and is sensitive to escaping 
particles such as high energy muons which leave only a minimum amount of energy in 
the calorimeter. 

• P^"^"-*- defined for events with at least one jet as the projection of the vector P^°'^° per- 
pendicularly to the jet axis. For events containing more than one jet the largest among 
all such projections is taken. In channels with missing neutrino signatures, a substantial 
pcaio± indicates that the missing momentum is not just due to fluctuations of the hadronic 
energy measurement. 



The selection criteria adapted to the different event topologies are described in subsec- 
tions ^yT| for excited fermions and|5]^ for excited quarks. For each of the possible decay chan- 
nels the number of selected events are compared to the Standard Model expectations. The 
errors given correspond to the statistical and the systematic uncertainties added in quadrature. 
A description of the systematic uncertainties can be found in section ^ 



5.1 Excited Leptons 



For e* and v* decays without muons or neutrinos in the final state, all particles are detected 
besides fragments of the proton. For these channels a cut 35 < [E — P^) < 65 GeV is applied 
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to reject photoproduction events where one jet is misidentified as an electron or a photon. The 
selection criteria for final states with muons are described in subsection |5.1.6| . For channels 
involving a neutrino P|,"'° > 20 GeV and {E — P^) < 50 GeV are required, with the additional 
cut P^^'-"-^ > 10 GeV when containing jets. For the channels with a Z or W boson in the 
final state which decays via W — > qq', Z qq or Z ee, a reconstructed invariant mass 
compatible with the boson mass within 20 GeV is imposed. This simple fixed size interval 
corresponds to about three times (seven times) the experimental mass resolution in the case of 
hadronic decays (leptonic decays). A similar mass cut is also imposed in the case of the decay 
chain u* eW ; W ^ euhy profiting from kinematic constraints (subsection p. 1.2p . 

5.1.1 The e* — > 67 channel 

The e* 67 decay mode is the key channel to search for e* because of its very clear signature 
and large branching ratio. The analysis starts from a sample of events with two electromagnetic 
clusters in the LAr calorimeter. The main sources of background are the QED Compton process, 
NC DIS with photon radiation or a high energy 7r° in a jet and the two-photon (77) production 
of electron pairs. Since about half of the cross-section is expected [0] in the elastic channel 
ep — > e*p, the analysis is split into two parts. The first (henceforward called "elastic") is 
dedicated to the search for e* produced elastically or quasi-elastically, the second ("inelastic") 
concerns the inelastic part of the e* cross-section ep — > e*X. 

• Elastic channel 

In this case the signature consists of only two electromagnetic clusters and no other sig- 
nals in the calorimeters. The clusters are required to have energies Ei> 5 GeV and angles 
6i < 150" {i = 1, 2), with a total energy sum Ei + E2 above 20 GeV. Because only high 
invariant masses of excited leptons are of interest, a cut M12 > lOGeV on the invariant 
mass M12 calculated from the two electromagnetic clusters is applied. The empty detec- 
tor condition consists in a cut Etot — Ei — E2 < 4 GeV, where Etot is the total energy 
deposited in the calorimeters. This cut strongly suppresses the NC DIS background. The 
remaining source of background is the elastic QED Compton scattering. 
After applying the above cuts 428 events remain. The expectation from Standard Model 
processes is 424 ± 10 (418.9 QED Compton and 5.1 77). 

• Inelastic channel 

Complementing the elastic analysis, we select here events with Etot — Ei — E2 > 4- GeV 
and require two high E^ electromagnetic clusters (Et^ > 20 and Et2 > 10 GeV and 
d-i < 150°). To reduce the NC DIS contribution with a high energy 7r° in a jet, a cut 
is applied on the multiplicity of tracks {ritracks < 2) in the direction close to that of the 
cluster with the second highest energy. 

With these criteria 150 candidates are selected, the expected Standard Model background 
is 158 ± 13 events (154.7 NC DIS events and 3.3 77). 

The combined efficiencies of the two selections are better than 80% for an e* mass above 75 
GeV. They are listed in table [1], as well as those of all other /* decay channels described later. 
The distributions of the measured "e — 7" invariant mass M12 are shown together with their 
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Table 1: Selection efficiencies in % for 
different decay modes of the excited lep- 
tons /* for masses Mi* ranging from 50 
to 250 GeV. The values given for the 
e* 67 correspond to the combined ef- 
ficiencies of the elastic and inelastic analy- 
ses. 



Ml. (GeV) 


50 


75 


100 


150 


200 


250 


e* 67 


73 


85 


87 


86 


83 


79 


v* 


24 


41 


47 


37 


32 


21 



Ml* (LreV) 


1 1 r> 


1 or\ 


1 <r> 
1 jU 


ZUU 




e eZ^ee 




/D 


/ / 


1 1 


/D 


e ^ eZ^^^ 






A Q 


J / 


ZO 


e* ^ eZ_^j.p 


< 1 


16 


72 


78 


80 




22 


44 


46 


40 


36 




27 


37 


44 


40 


34 




11 


29 


40 


34 


29 


Z/* ^ vZ^qq 


11 


25 


41 


46 


37 


V* ^ el^_^e. 


47 


45 


42 


38 


32 




16 


27 


40 


35 


35 




38 


39 


38 


34 


28 



a) Search for e — > e y : elastic analysis 

■ HI data 

□ QED Compton 




50 100 150 200 
e-y invariant mass (GeV) 



^^10^ 

o 

* 

^10^ 



§10 



10 
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b) Search for e — > e y : inelastic analysis 

■ HI data 

NC DIS 




50 100 150 200 
e-y invariant mass (GeV) 



Figure 2: Invariant mass spectrum for (a) the elastic and (Z?) the inelastic e* — > 67 analyses. The 
invariant mass is calculated for an event by combining the four-momenta reconstructed from the 
two electromagnetic clusters. Square symbols correspond to the data and the histograms to the 
expectations from the different Standard Model processes. 
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corresponding expected background in fig. ^ and fig. ^ for the elastic and inelastic analyses 
respectively. The expected mass resolution {a of a Gaussian fit) for the reconstruction of an e* 
mass of Afg. = 150 GeV is 4.2 GeV. The number of events with M12 above 50 GeV is 53 in 
the elastic analysis and 68 in the inelastic analysis, in good agreement with the predictions from 
the Standard Model of 52.5 ± 1.5 and 77.4 ± 8.2 respectively. 

5.1.2 The e* — > eZ and i^* eW, vZ channels with Z — > ee and W ^ e.v 

The selection for these channels starts from a sample of events with two high Ej- electromag- 
netic clusters {Er^ > 20 and Et^ > 10 GeV and 6, < 150°). 

• Events with three electromagnetic clusters 

In the case of the e* eZ^ee channel, a third electromagnetic cluster with Et.^ > 5 GeV 
and 63 < 150° is required. Furthermore at least one pair of the three electromagnetic 
clusters should have an invariant mass compatible with the Z mass. To reduce NC DIS 
background, events having a jet with E-r^^^ above 15 GeV are rejected. After this selection 
1 event is left while 0.9 ± 0.4 background events are expected (0.1 from NC DIS and 0.8 
from 77 processes). 

• Events with two electromagnetic clusters and missing Pt 

In these cases, the four-momentum of the escaping neutrino is deduced by imposing 
transverse momentum conservation and the (E — Pz) constraint. In the case of the 
u* — ^ eW^eu and v* i^Z^^.^. decay channels, events are selected with the invariant 
mass from the u and electromagnetic cluster (for W tagging) and from the two electro- 
magnetic clusters (for Z tagging) to be compatible with the corresponding boson mass. 
No candidate is found. The expectation from Standard Model processes is 0.25 ±0.11 
events (0.05 77, 0.1 NC DIS, 0.1 W ^ eu) for the u* eW channel, and 0.020 ± 0.005 
W ^ ev events for the u* vZ channel. 

5.1.3 The e* eZ, vW and v* vZ, eW channels with Z,W ^ qq' 

The analysis for these channels uses a subsample of events with at least two jets each having a 
transverse energy greater than 20 GeV and a polar angle greater than 10°. The jet-jet invariant 
mass must be compatible with the relevant boson mass. When more than two jets are found 
in an event, the pair of jets which has an invariant mass closest to the relevant boson mass is 
selected. This subsample is dominated by photoproduction and NC DIS events. 

• Events with two high Et jets and one electron 

The channels e* eZ^gg and u* eW^qq are characterized by two high Et jets and an 
electron. Background events are expected from NC DIS. Candidates are selected if they 
have an electron with Pt^ > 15 GeV and with 10° < 9^ < 90°. This cut on the transverse 
momentum of the electron induces an efficiency loss towards low e* masses, already 
sizeable close to the Z mass (see table |1]). The cut on the polar angle of the electron 
discriminates the signal, where the lepton is mainly emitted in the forward direction due 
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the high /* mass, from NC DIS background where the electron is mainly scattered in the 
backward region. 18 events are found in the search for the channel v* eW with an 
estimated background of 17.2 ± 4.8 events (16.6 NC DIS, 0.4 7p and 0A6W qq'). 14 
events survive the selection criteria for the decay e* eZ, in comparison to a background 
of 12.3 ± 3.4 events (12 NC DIS, 0.14 7p and 0.16 W qq'). Fig. gshows the invariant 
masses M^jj of the two jets and the electron. For a /* mass of 150 GeV, the expected mass 
resolutions are 8.6 and 13.8 GeV for the e* eZ and eW channels, respectively. 
The u* mass resolution is worse than the e* one, because the recoil jet is in some cases 
wrongly taken as one of the jets associated to the W decay. No excess of events is found 
compared to the Standard Model expectation. 



a) Search for v — > e W 



■ HI data 

m NCDIS 




50 100 150 200 

e-.jet-jet invariant mass (GeV) 



b) Search for e — > e Z 



■ HI data 

NC DIS 




100 150 200 

e-jet-jet invariant mass (GeV) 



Figure 3: Invariant mass spectrum for the (a) u* — > eW^qq and (b) e* — > eZ^gg searches. The 
invariant mass is calculated for an event by combining the four-momenta reconstructed from 
the electromagnetic cluster and the two jets. Square symbols correspond to the data and the 
histogram to the expectation from the NC DIS process. 



• Events with two high Et jets and missing Pt 

In the searches e* uW^qg and u* — > vZ^gg, the main background which is due to 
CC DIS interactions is suppressed by the P^"'°, {E - P,) and P^^°^ cuts. The NC DIS 
background is reduced by rejecting events possessing an electromagnetic cluster with an 
energy above 5 GeV. Three events survive. The background expectation is 3.3 ± 0.6 
events (3 CC DIS, 0.3 7p) for the e* uW^gg channel and 2.1 ± 0.8 events (2 CC and 
0.1 7p) for the v* — > vZ^gg channel. 
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5.1.4 The u* vj channel 

For this analysis events containing one photon with Et^ > 20 GeV and 10° < 9^ < 90° and 
satisfying the P|"'" and {E — Pz) cuts are selected. As in subsection |5.1.3| , the aim of the 



9 cut {9^ < 90°) is to further reduce the NC DIS background for which the electromagnetic 
cluster is predominantly reconstructed in the backward region. The final state for the signal 
contains also in most of the cases a recoil jet, due to the u* production through a t— channel 
W boson exchange. Hence the final selection criteria are one jet with Et^^^. > 5 GeV and no 
electron found with an energy above 5 GeV. No candidate is left. The expected Standard Model 
contribution is 1.0 ± 0.7 events and is dominated by CC DIS events. 



5.1.5 The e* — > eZ^i^p channel 

This channel is characterized by one high electron and missing Pt in the detector. For 
this channel the non ep background (cosmic rays and halo muons) is severe and hence the 
minimum P^"'° requirement is increased from 20 to 25 GeV. A requirement of large transverse 
momentum for the electron is also necessary to reduce the background, so events with a high 
Pt electron (Pt^ > 20 GeV with 10° < 9^ < 150°) are selected. Events with either another 
electromagnetic cluster of energy greater than 5 GeV or with a jet with Et above 15 GeV are 
rejected. This selection finds 1 event for an expected background of 2.7 ± 0.4 events (1.3 NC 
DIS, 0.7 W eu, 0.6 CC DIS and 0.1 77). 



5.1.6 The u* eW^^^ and e* eZ^^^ channels 

The search in these channels starts from a subsample of events including at least one muon 
candidate found at a polar angle greater then 10° with a transverse momentum above 10 GeV. 

• Events with one muon and an electron 

For the v* eW^^y analysis we require P™'° > 25 GeV and a high P^ electron (P^^ > 
20 GeV). No events are left after this selection and the total expected background which 
survives these two cuts is 0.31 ± 0.05 events dominated by the muon pair production in 
77 interactions (0.28 events) with a small contribution (0.03 events) from single W 
production. It has been checked that applying the more stringent cuts of P^ ] reduces the 
77 background by a factor of ~10. 

• Events with two muons and an electron 

The signature of the e* — > eZ^^^ channel consists of two muons plus an electron. Here, 
contrary to the preceding one muon case, no neutrino is expected. However as the energy 
deposited in the calorimeter by the two muons is small, a cut P™'° > 15 GeV is applied. 
With the requirement of two identified muons no events are left. The background due to 
the 77 nn process, is 0.35 ± 0.05 events. 
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5.2 Excited Quarks 



In a q* production process at HERA one would expect that the scattered positron be often unseen 
in the detector since the process is dominated by 7 exchange at small values of Q^. Hence no 
restriction is imposed on the value of (E — P^). The selection criteria for channels with one or 
two muons are described in subsection |5.2.5| . For channels with a neutrino, the cuts P™'° > 20 
GeV and P|"'°-L > 10 GeV are applied. A compatibility within 20 GeV of the Z or W boson 
mass is imposed here also in the case ofW^ qq', Z qq or Z ^ ee decays. 



5.2.1 The q* q'y channel 

The final state for this channel is characterized by one photon and one jet. The photon require- 
ments are those described in subsection |5.1.4] : one photon with {Et^ > 20 GeV, 10° < 9^ < 
90°). In addition we require one jet with -^Tj^t > 15 GeV and Ojet > 10° and no electron 
with an energy above 5 GeV in the LAr calorimeter. The main background sources are pho- 
toproduction processes with prompt 7 production or events with high energy 7r°, and NC DIS 
events if the track of the scattered electron has not been reconstructed. 35 events are selected 
compared to a background estimation of 36 ± 5 events (2.5 NC DIS and 33.5 7p). Fig. ^ shows 
the distributions of the invariant 7-jet masses for the data and the expected background. The 
expected invariant mass resolution for a q* mass of 150 GeV is 6.6 GeV. The two events with 
an invariant mass (7-jet) above 150 GeV are very likely NC DIS events, both just surviving the 
criteria on the quality of the nearest track or on its distance to the electromagnetic cluster. The 
efficiency of this selection is listed in table 0, as well as those of all q* decay channels described 
later. 



Figure 4: Invariant mass 
spectrum for the q* — > 
57 search. The invariant 
mass is calculated for an 
event by combining the four- 
momenta of the photon and 
the jet. Square symbols 
correspond to the data and 
the histograms to the expec- 
tations from different Stan- 
dard Model processes. 
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Table 2: Selection efficiencies in % 
for different decay modes of the excited 
quarks q* for masses Mg. ranging from 50 
to 250 GeV. 



Mg* (GeV) 


50 


75 


100 


150 


200 


250 


q* g7 


22 


33 


36 


42 


41 


40 



Mg, (GeV) 


110 


120 


150 


200 


250 


q* qZ^ee 


16 


35 


40 


47 


41 


q* qZ^^^ 


7 


16 


25 


32 


32 


q* qZ^^q 


< 1 


3 


31 


47 


44 


q* qW^eu 


31 


38 


42 


44 


38 


q* qW^f,^ 


6 


24 


41 


41 


37 


q* qW^qq 


5 


16 


42 


51 


47 



5.2.2 The q* qZ^^^ channel 

Starting from the subsample of events with at least two electromagnetic clusters described in 
subsection |5.1.2| , it is required in addition that there be a jet with Et^^^ > 15 GeV and 9 jet > 10° 



and that the invariant mass calculated from the two electromagnetic clusters be compatible 
with the Z mass. No events survive these criteria, a fact consistent with the Standard Model 
expectation of 0.65 ± 0.53 NC DIS events. 



5.2.3 The q* qW^eu channel 

For this channel events retained by the P^°-^", (E — P^) and P^°-^°-^ cuts and with only one high 
Pt electron {Pt, > 15 GeV and 10° < Oe < 150°) and a jet (Et^^^ > 15 GeV and Ojet > 10°) 
are selected. The main sources of background in this channel are NC and CC DIS processes 
and W production. One event survives this selection and the number of events expected from 
the Standard Model is 1.10 ± 0.35 events, equally shared between the NC DIS and W eu 
processes. 



5.2.4 The q* qW-^qq> and q* qZ^qq channels 

The final state in these channels contains three high Et jets and the main backgrounds are 
photoproduction and NC DIS processes. We require three jets with a polar angle above 10° 
and transverse energies greater than 30, 25 and 15 GeV, respectively. Furthermore, in these q* 
decays, the jet with lowest transverse energy often originates from the boson (W/Z) decay, when 
the q* mass is above 150 GeV. Events are kept only when the jet-jet invariant mass calculation 
which is the nearest to the W/Z mass includes this jet. For the channel q* — > qZ^gg, the loss 
of efficiency due to this requirement varies between 30% for a q* mass of 150 GeV to less than 
one percent when the q* mass is equal to 250 GeV. In the q* qW^qq case the losses are 
smaller. We require that this invariant jet-jet mass be compatible with the boson mass. These 
cuts select 39 and 32 events for the q* qW^qq and the q* qZ^qq channels, respectively. 
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These numbers are to be compared to the Standard Model expectations of 45.3 ± 17.3 events 
(30.4 7p, 13.2 NC DIS and 1.7 W) for the q* qW analysis, and of 25.3 ± 9.1 events (17.6 
7P, 6.4 NC DIS and 1.3 W) for q* qZ. The NC DIS expectation has been calculated using 
the generator RAPGAP, in which leading log parton showers are used to model QCD radiations. 
The DJANGO generator which uses the CDM to simulate QCD effects is not able to describe 
the data in this particular phase space domain. A similar observation in the measurement of 
2-jet rates in DIS has been reported previously [|T7[]. Detailed investigations of discrepancies 
between different QCD cascade models are underway []38|]. The shapes of the invariant three-jet 
mass distributions are in good agreement with the Standard Model expectations, as can be seen 
in fig. |[ For a q* mass of 150 GeV, the expected resolution on the 3-jet invariant mass is 9 
GeV. 



a) Search for q — > q W 



■ HI data 
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b) Search for q q Z 

■ HI data 



■ w 

□ yp 
H NC DIS: 



50 100 150 200 250 

3-jets invariant mass GeV 



Figure 5: Invariant mass spectrum for (a) the q* — * qW^qq and {b) the q* qZ^qq analyses. 
The invariant mass is calculated for an event by combining the four-momenta of the three jets. 
Square symbols correspond to the data and the histograms to the expectations from different 
Standard Model processes. 



5.2.5 The q* qW^n.i, and q* qZ^^^^^ channels 



These analyses use the muon subsample described in subsection |5.1.6| , together with the cuts 



pcaio ^ 25 GeV and P^'^'-" > 15 GeV for the one muon and two muon searches, respectively. 

• Events with one muon and a jet 

In the q* qW^^j,^ search, events including a high Et jet {ET^^t > 25 GeV and 9jet > 
10°) are selected. To reduce the 77 — > fi^ background, events with more than one isolated 
muon are rejected and an acoplanarity cut A0(/i — jet) < 175° in the transverse plane is 
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applied. Three events are found. They correspond to the events labelled muon-2, muon-4 
and muon-5 observed already in a previous search for events with high Px leptons and a 
large missing transverse momentum [RH|. The background (0.41 ± 0.03 events) is here 



due to W production (0.31 events) and 77 muon pair production (0.1 events). 
Events with two muons and a jet 

The signature of the q* qZ^^^ channel consists of two muons plus a jet. When request- 
ing two muons and applying a P^"'" > 15 GeV cut no events remain. The background is 
equal to 0.35 ± 0.05 events from the 77 nn process. 



6 Limits on Excited Fermion Production 



6.1 Upper Limits on Cross-sections 



A summary of the number of events surviving the selection cuts in the various channels is 
given in table |3| for the excited lepton decay channels and in table ^ for the excited quark decay 
channels. The uncertainties taken into account on the background determination are listed in 
the following. 



• The statistical error of the Monte Carlo generations. 

• An uncertainty on the absolute electromagnetic energy scale ranging from ±0.7% in the 
central part of the LAr calorimeter to ±3% in the forward part. 

• An uncertainty of ±4% on the energy of the jets due to the uncertainty on the calibration 
of the calorimeters for hadronic showers. 

• A 7% uncertainty on the predicted DIS cross-sections coming mainly from the lack of 
knowledge on the proton structure (see detailed discussion in [pO[]). 

• An uncertainty of ±10% on the expectation for the 2-jet cross-section estimated by com- 
paring leading order and next-to-leading order Monte Carlo simulations. In the same way 
an uncertainty of ±15% on the three-jet cross-section was determined by a comparison 
to data of either a Monte Carlo with 0{as) QCD matrix elements which approximates 
the higher order emission of partons using the concept of parton showers, or perturbative 
QCD calculations to order O(a^) which produces an exact leading order calculation of 
the three parton final state pT]]. 

• An uncertainty on the estimation of the radiative CC DIS background for the u* — > wj 
channel, coming from the fact that photon radiation from the quark lines is not fully taken 
into account in our CC DIS simulation. Calculations show that the negative interfer- 
ence term between photon radiation from the electron and quark lines could decrease the 
radiative CC DIS cross-section by an amount of 70%. 

• An overall systematic error of 1.5% on the luminosity. 
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The observed number of events are compared to the expected Standard Model background 
in tables |3| and ^. Good agreement is found for all channels, except for a slight excess, already 
quoted in reference [^, observed in the q* qW^^^ channel where 3 events are observed 
for a Standard Model expectation of 0.41 ± 0.03 events. However, combining the contributions 
from the three W decay channels, no significant deviation to the Standard Model prediction 
remains within the present analysis. 



Channel 


selection criteria 


events 


11 1 

background 


e* — > 67 (cl. an.) 


2 em clusters only 


J J 


jZ.j ± 1 .J 


e* 67 (inel. an.) 


2 high Et em clusters 


68 


77.4 ± 8.2 


e* eZ^ee 


3 em clusters 


1 


0.9 ± 0.4 


e* eZ^f,^ 


2 muons + 1 electron 





0.35 ± 0.05 


e* eZ^yy 


1 electron + P^''' 


1 


2.7 ± 0.4 


e* eZ^gg- 


2 jets + 1 electron 


14 


12.3 ± 3.4 


e* vW^qq 


2jets + P^*^" 


3 


3.3 ±0.6 


V* — > z/7 


1 photon + P^''' 





1.0 ±0.7 


V* vZ^ee 


2 em clusters + P™"** 





0.020 ± 0.005 


V* uZ^qq 


2jets + Pf^*^" 


3 


2.1 ±0.8 


V* eW^eu 


2 em clusters + P™'''' 





0.25 ±0.11 




1 muon + 1 electron 





0.31 ±0.05 


U* ^ eW^qq 


2 jets + 1 electron 


18 


17.2 ±4.8 



Table 3: Number of events observed in the various e* and u* decay channels and the corre- 
sponding Standard Model expectation and total uncertainty on the mean expectation. It should 
be noted that these numbers correspond to different invariant mass intervals, as the effective 
mass threshold depends on the channel. 



No evidence was seen for either excited leptons or quarks in any of the channels. Therefore, 
upper limits on the product of the f * production cross- section and the decay branching fraction 
have been derived. These limits are determined at a Confidence Level (CL) of 95% as a function 
of the excited fermion mass. A mass window is shifted over the whole mass range in steps of 
5 GeV. The width of each window is chosen according to the resolution for the corresponding 
mass. When combining several decay channels, for each decay channel k and in each mass 
interval, the number of observed events Uk, the number of expected background events bk and 
efc, the product of the efficiency times branching ratio of the channel, are calculated and used to 
determine the value A of the upper limit for the signal such that: 

/•A roo ^ 1 

CL= p{a)da/ p{a)da ; p{a) = 11— {eta + bkT' exp-^''^''+''^^ (4) 

Jo Jo Y ^k' 

where a is the Poisson parameter of the signal. For a single decay channel this is identical to 
the Bayesian prescription given by the Particle Data Group [^. For the background estimation 
and the selection efficiency, statistical and systematic errors are taken into account by folding 
Gaussian distributions into the integration of the Poisson law used to determine the limit. 

Because the branching ratios of the hadronic decay of the W or Z bosons are dominant, 
the limit for the f * decaying into a fermion and a W or a Z mainly depends on final states 
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Channel 


selection criteria 


events 


background 


q* g7 


1 photon + 1 jet 


35 


36 ±5 


q* qZ^ee 


2 em clusters + 1 jet 





0.65 ±0.53 


q* qZ^^^ 


2 muons + 1 jet 





0.35 ± 0.05 


q* qZ^qq 


3 jets 


32 


25.3 ±9.1 


q* qW^ev 


1 electron + 1 jet + P™*'*'' 


1 


1.10 ±0.35 


q* qW^i^y 


1 muon + 1 jet 


3 


0.41 ± 0.03 


g* ^ qW^,, 


3 jets 


39 


45.3 ± 17.3 



Table 4: Number of events observed in the various q* decay channels and the corresponding 
Standard Model expectation and total uncertainty on the mean expectation. It should be noted 
that these numbers correspond to different invariant mass intervals, as the effective mass thresh- 
old depends on the channel. 



with at least two jets. So the error is dominated by the uncertainty on the absolute calibration 
of the calorimeters for hadronic clusters and the uncertainty on the expected 2-jet and 3-jet 
cross-sections. 

The limits on the product of the f * production cross-section and the decay branching frac- 
tion are shown in fig. ^ ^ and ||. In all three cases the lowest limits are obtained using the 
electromagnetic decay channels. For the derivation of these limits it is assumed that the natu- 
ral width of the f * is much smaller than the experimental mass resolution. For masses above 
120 GeV, the values of these limits are below 0.2 pb for the e*, 0.5 pb for the u* and 0.4 pb for 
the q* productions. These results improve by an order of magnitude earlier HI results [|^]. 
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Figure 6: Upper limits at 95% Con- 
fidence Level on the product of the 
production cross-section a and the de- 
cay branching fraction BR for excited 
electron e* in the various electroweak 
decay channels, 67 (full line), eZ 
(dashed line) and vW (dotted-dashed 
line) as function of the excited elec- 
tron mass. The different decay chan- 
nels of the W and Z gauge bosons are 
combined. Areas above the curves are 
excluded. 
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Figure 7: Upper limits at 95% Con- 
fidence Level on the product of the 

production cross-section a and the de- 
cay branching fraction BR for excited 
neutrino v* in the various electroweak 
decay channels, 1/7 (full line), vZ 
(dashed line) and eW (dotted-dashed 
line) as function of the excited neu- 
trino mass. The different decay chan- 
nels of the W and Z gauge bosons are 
combined. Areas above the curves are 
excluded. 
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Figure 8: Upper limits at 95% Con- 
fidence Level on the product of the 
cross-section a and the electroweak 
decay channel branching ratio BR for 
excited quark production in the elec- 
tromagnetic (full line), the Z (dashed 
line) and W (dotted-dashed line) de- 
cay channels as function of the excited 
quark mass. The different decay chan- 
nels of the W and Z gauge bosons are 
combined. Areas above the curves are 
excluded. 
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6.2 Upper Limits on Coupling Parameters 

Assuming fixed numerical relations between / and /', the cross-section depends only on //A 
and Mf*, and thus constraints on //A can be derived. Conventional assumptions are / = /' or 
/ = — /'. From the coupling constant relations (see equation 2) it can be seen that the coupling 



20 



of the e* to the 67 mode is proportional to (/ + /'), and that of the u* to the 2/7 decay channel 
is proportional to (/ — /'). 

In fig. 1^, limits on the ratio //A are given for the e* for the hypothesis / = /'. We do 
not consider the case / = — /', because the coupling constant c^e*e would be equal to and 
the production cross-section of the e* is then very small. Upper limits for //A ranging from 
7 X 10-^ GeV-i to 10-2 GeV^^ are obtained at 95 % CL for an e* mass ranging from 50 GeV 
to 250 GeV. 

In fig. [T^, limits on the ratio //A are given for v* , assuming / = — /' and / = /'. These 
two assumptions correspond to very different u* branching ratios, as shown in table |] for the 
example of two u* masses. In particular when f = f the u* 67 has a branching ratio equal 
to 0. Somewhat better limits are obtained when / = — /'. The values of the limits for //A vary 
between 3 x 10^'^ to 10^^ GeV^^ for an mass ranging from 50 to 200 GeV. These upper 
limits are conservative for masses above ~ 170 GeV where the narrow-width approximation 
underestimates the total u* production cross-section. 

Assuming //A = 1 /Mi*, masses below 223 and 1 14 GeV are excluded at 95 % CL for the 
e* (/ = /') and u* (/ = — /') production, respectively. 
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Table 5: Branching ratios in % of the 
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Limits for the q* assuming f = f and only electroweak couplings (i.e. fs = 0) are shown 
in fig. |rT| The exclusion limits for q* masses between 50 to 250 GeV corresponds to values 
of //A between 9 x IQ-^ and 2 x lO^^ GeV"\ Assuming //A = l/M^., masses below 188 
GeV are excluded at 95 % CL. The /s = assumption allows to make a study which is comple- 
mentary to the analysis done by the CDF experiment M at the Tevatron. At a pp collider 



excited quarks are produced in a quark-gluon fusion mechanism which requires fs 7^ 0. As- 
suming A = Mg*, CDF excludes excited quarks in the mass range 80-300 GeV for f = f = fs 
values greater than 0.2 and up to 760 GeV if / = /' = = 1. 

The complementarity of the HI results to those of CDF is illustrated on fig. where a compar- 
ison of the exclusion domains of / = /' values, obtained by CDF and HI, is shown for different 
hypotheses on the fs value. As soon as fs is smaller than ~ 0.1 and for Mg* < 130 GeV, our 
analysis probes a domain not excluded by Tevatron experiments. 

LEP experiments have also reported searches for excited leptons at e+e^ center of mass 
energies up to ^/s = 189 GeV ^ The results for excited leptons produced in pairs. 
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Figure 9: Exclusion limits on cou- 
pling constants at 95% Confidence 
Level as a function of the mass of 
the excited electron. The assumption 
/ = /' is made for the different de- 
cay channels (full, dashed and dotted- 
dashed lines) and for all decay chan- 
nels combined (dotted line). Values 
of the couplings above the curves are 
excluded. The light grey area corre- 
sponds to the exclusion domain ob- 
tained by HI in this analysis. The dark 
grey area is excluded by the L3 exper- 



ts 
O 



^0 



10 



iment [45]. 



> 
O 



a)f=-f' 

— v*-^ V X 

v^-^vZ 

v*-^ eW 

V all channels 



10 



10 



10 




100 150 

V* Mass 



200 

(GeV) 



10 



10 





1 i 1 1 1 1 1 1 i 


1 1 1 1 1 1 1 . 

HI 


— 






_ 




f = f ' ■ 


e. 






e 


all channels 






1 , . . . 1 . . , 


, 1 .... 1 . 



-4 

50 



100 



150 200 

e Mass 



250 

(GeV) 



b)f = f' 



V -> V Z 
v*-> e W 

V ^ all channels 



10 



10 



10 



HI 




Figure 10: Exclusion limits on coupling constants at 95% Confidence Level as a function of the 
mass of the excited neutrino. The assumptions f = —f and f = f are made for figures (a) and 
(b) respectively. The results for the different decay channels are shown separately (full, dashed 
and dotted-dashed lines) and for all decay channels combined (dotted line). Values of the cou- 
plings above the curves are excluded. The light grey area corresponds to the exclusion domain 



obtained by HI in this analysis. The dark grey area is excluded by the L3 [4-5] experiment. 



when the limits are independent of / and /', as well as results for single e* and u* production 
assuming / = /' or / = — /' are included in fig. ^ and [1^. In [ p7| ] results independent from 
a hypothesis on the relation between / and /' are also given. Fig. [TT] shows the result of a q* 
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Figure 11: Exclusion limits on cou- 
pling constants at 95% Confidence 
Level as a function of the mass of 
the excited quark, assuming f = f 
and fs = 0. The results for the 
different electroweak decay channels 
are shown separately (full, dashed 
line and dotted-dashed lines) and for 
all decay channels combined (dotted 
line). Values of the couplings above 
the curves are excluded. The light 
grey area corresponds to the exclu- 
sion domain obtained by the HI ex- 
periment in this analysis. The dark 
grey area is excluded by the DEL- 
PHI experiment [ |4^ ] assuming that 
the branching ratio of the q* q'j is 
equal to 1 . 
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Figure 12: Exclusion limits on / 
values at 95% Confidence Level as a 
function of the mass of the excited 
quark, assuming A = M{q*), f = f 
and for different fg values. Exclusion 
limits from CDF (the 2 right curves) 
for 2 fs values have been derived from 
table 1 of reference [M. Values of 



the couplings above the curves are ex- 
cluded. 



A = M(q') f = f 

HI 




search [ ]4^ ] at LEP, assuming a branching ratio of 1 for q* 
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7 Summary 



Using e'^p data taken from 1994 to 1997 corresponding to an integrated luminosity of 37 pb^ , a 
search for production of excited fermions has been carried out. No evidence has been found for 
excited electrons, excited neutrinos or excited quarks for decays into any of the gauge bosons 
'y, Z,W and Standard Model fermions. 

New limits for the production of excited fermions have been obtained, which improve previ- 
ous HI results by a factor 10 and previous published ZEUS results [Q| based on 4 times smaller 
integrated luminosity. For masses above 180 GeV, i.e. in a domain extending beyond the kine- 
matic reach of LEP, compositeness scales of //A in the range of 1.5 x 10^'^ to 2 x 10^^ GeV^^ 
are excluded from the search for excited electrons, and 5.6 x 10^^ to 0.32 GeV"^ for excited 
neutrinos. Assuming //A = 1/Mf*, excited fermions with masses below 223, 114, and 188 
GeV, for e*, u* and q* productions, respectively, are excluded. The results obtained on q* pro- 
duction via electroweak couplings are complementary to the results obtained at the Tevatron pp 
collider where q* production via strong coupling is investigated. 
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